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ABSTRACT 

Through  use  of  the  low -density,  hypervelocity  wind  tunnel  of  the  von 
Karman  Gas  Dynamics  Facility,  drag  of  spheres  has  been  measured  under 
hypersonic,  cold- wall  (Tw  «  T05,  support-free  conditions  in  a  nunreaeting 
(vibration  frozen)  flow.  Data  were  obtained  for  a  nominal  free -stream 
Mach  number  of  11  and  for  Reynolds  numbers  from  1  to  10  based  on  condi¬ 
tions  immediately  downstream  of  the  bow  shock  and  sphere  diameter. 
These  data  were  supplemented  by  measurements  at  a  nominal  Mach  num¬ 
ber  of  10  where  a  conventional  balance  was  used,  and  Reynolds  numbers 
downstream  of  the  shock  as  high  as  104  were  investigated  in  the  cold -wall 
condition. 

The  experimental  results  have  been  analyzed  both  from  the  point  of 
view  of  continuum  flow  with  second- order  viscous  effects,  and  from  the 
standpoint  of  a  noncontinuum  concept,  taking  account  of  first  collisions 
between  reemitted  and  free-stream  molecules.  In  the  former  case,  a 
drag  coefficient  representation  of  the  form 

Cn  “  Cpjirfiacid  +  +  K,/Ret 

is  shown  to  fit  the  data  very  closely  when  K,  and  K.,  are  allowed  to  be 
functions  of  free-stream  Mach  number  and  wall -to-total  enthalpy  ratio. 
Data  from  investigations  at  the  University  of  California,  University  of 
Toronto,  and  Jet  Propulsion  Laboratory  have  been  used  to  enlarge  the 
present  study  and,  in  particular,  to  support,  the  evaluation  of  K,  which 
represents  the  influences  of  vorticity,  curvature,  thick  boundary  layer, 
slip,  and  temperature  jump.  In  all  cases,  K2  was  found  to  be  negative 
in  sign. 

The  new  first-collision  analysis  is  numerically  indeterminant  be¬ 
cause  of  the  lack  of  a  method  for  explicit  calculation  of  mean  free  path 
in  real,  polyatomic  gases  with  consideration  of  intermolecular  forces. 
However,  the  form  of  the  derived  equation  for  (JD/Cn,ree  ooiecai,,  appears 
to  fit  the  experimental  data  over  a  considei’able  range  of  Knudsen  num¬ 
bers  if  a  free  constant  is  used  in  the  expression  for  mean  free  path.  It 
is  considered  that  this  constant  is  influenced  by  the  type  of  intermolecu¬ 
lar  and  gas-solid  interactions  occurring.  Fully  accommodated,  diffusely 
reflected  surface  interaction  was  assumed  in  the  analysis. 
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NOMENCLATURE 


A 

a 

B 

b 

Clt  Clf  C, 


Cn 

cn(ro 

Cni 

Cp 

Cv 

D 

a 

dA 

do> 

-Ei  (-x, 
erf  x 
G  ( v/  v  ) 

8 

Hw 

h 

K 


Ki.K, 


Kn 

k 

i 


M 


m 

N(i) 


v(r/Aw)2'-  (r/AwT* 

Speed  of  sound,  \/y/?T 

r/Aw  -  R/Aw 

Constant 

Constants 

Drag  coefficient,  defined  by  Eq.  (3) 

Free-molecule  drag  coefficient,  defined  by  Eq.  (53) 

Inviscid  drag  coefficient 

Specific  heat  of  gas  at  constant  pressure 

Specific  heat  of  gas  at  constant  volume 

Drag  of  sphere 

Diameter  of  sphere 

Elemental  area  of  surface 

Incremental  solid  angle 

Generalized  exponential  integral,  f.  dx 

Error  function 

Function,  defined  by  Eq.  (43) 

Acceleration  due  to  gravity 

Enthalpy  parameter,  defined  by  Eq.  (7) 

Enthalpy 

Thermal  conductivity 
Constants 
Knudsen  number 
Boltzmann  constant 

Coordinate  measured  from  a  point  on  a  surface  to 
a  point  in  space 

Mach  number 
Mass  of  a  molecule 
Number  of  molecules 
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n 

n  (v) 


P 

Pr 

ft 

R 

Re 

r 

S 

T 

U 


X 

y 

a  ( v ) 


y 


Ax 


n 


e 

x 

A(  v) 
Aw 

A„ 


Number  density  of  molecules 

Number  density  of  molecules  in  velocity, 
position  coordinates 

A  point  in  space 

Prandtl  number 

Radius  of  sphere 

Gas  constant,  k/as 

Reynolds  number,  Vdp/fi 

Coordinate  from  center  of  sphere 

Molecular  speed  ratio,  U/ 2R T 

Temperature 

Velocity 

Velocity  of  a  molecule 

Mean  molecular  speed,  y/~BRT/n 

Mean  velocity  of  molecule  leaving  a  surface 

Axial  coordinate 

Vertical  coordinate 

Collision  probability  of  a  molecule  with  speed 
v ,  1/  A  ( v ) 

Ratio  of  specific  heats,  Cp/Cv 
Difference  between  axial  coordinates 

Collision  rate  density  between  wall  and  free- 
stream  moiecuies 

Angle  between  normal  to  dA  and  the  free- stream 
velocity  vector  (see  Fig.  10) 

Moan  frog  path 

Mean  free  path  of  molecules  with  speed  v 

Mean  free  path  of  a  molecule  reemitted  with 
speed  vw  relative  to  free-stream  molecules 

Mean  free  path  of  reemitted  molecules  relative 
to  free-stream  molecules 

Viscosity 


x 


Mass  density 

Effective  molecular  diameter 
Shear  stress 
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Angle  between  the  normal  to  <M  and  the  incident  or 
reemitted  path  (or  velocity  vector)  of  a  molecule 
(see  Fig.  10) 


SUBSCRIPTS 

o  Total  or  isentropic  stagnation  condition 

2  Condition  immediately  downstream  of  normal  shock 

i  Incident  molecules 

a  Scattered  molecules 

w  Wall  condition 

~  Free- stream  condition 
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1.0  INTRODUCTION 


It  is  well  known  that  the  aerodynamic  drag  of  spherical  bodies  in 
hypersonic  flow  is  more  than  doubled  when  the  flow  passes  from  a  con¬ 
dition  characterized  by  '’thin'1  boundary  layers  to  a  condition  of  free- 
molecular  flow.  The  first  condition  usually  is  identified  with  large 
Reynolds  numbers,  Re,  whereas  the  flow  parameter  associated  with  the 
second  condition  is  Knudsen  number,  Kn.  A  relation  between  these  param¬ 
eters  and  the  Mach  number,  M,  was  derived  by  von  Karman  who  found 

Kn  -  M/Re  (1) 

The  constant  of  proportionality  is  approximately  1.  49  for  air  at  moderate 
temperature  and  pressure. 

The  drag  coefficients  of  hypersonic  spheres  at  both  very  high  Knudsen 
numbers  (>  10)  and  very  low  Knudsen  numbers  (<  10~4)  are  known  to  at 
least  a  moderate  degree  of  accuracy,  although  there  is  uncertainty  regard¬ 
ing  the  accommodation  coefficients  in  free-molecular  flows.  In  the  inter¬ 
mediate  range  of  Knudsen  numbers,  where  neither  analysis  nor  experiment 
has  fully  answered  the  question,  significant  contributions  have  been  made 
by  both  approaches.  Some  of  the  more  recent  work  is  reported  in 
Refs.  1-16,  where  sources  of  additional  material  also  are  given. 

With  regard  to  the  high  Knudsen  number  or  near-free-molecular  flow 
regime,  the  obstacles  to  analytical  solution  arise  chiefly  from  the  lack 
of  exact  knowledge  of  the  behavior  of  molecules  after  impact  on  the  sphere 
and  the  mathematical  complexity  of  the  flow  mode1  when  interactions 
between  free-stream  and  reflected  molecules  become  important.  In  the 
continuum,  viscous-interaction  regime  of  small  Knudsen  numbers  the 
usual  difficulty  of  dealing  with  flows  characterized  by  thick  boundary  layers 
in  the  presence  of  vorticity,  curvature,  displacement,  slip,  and  tempera¬ 
ture  jump  presents  the  main  obstacle. 


The  reason  that  experimental  data  have  not  been  presented  for  the 
entire  intermediate  regime  is  basically  the  lack  of  facilities  capable  of 
duplicating  all  the  conditions  of  flow  desired  for  the  investigation.  None¬ 
theless,  experimental  data  provide  the  more  significant  part  of  the  know¬ 
ledge  of  sphere  drag  in  the  region  10"“  <  KnM  <  10  ,  where  Kn^  is  based  on 
free-stream  conditions  and  sphere  diameter. 


Manuscript  released  by  authors  October  1962. 
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The  Low- Density,  Hypervelocity  (LDH)  Wind  Tunnel  (Ref.  7)  of 
the  von  Karman  Gas  Dynamics  Facility  (VKF),  Arnold  Engineering 
Development  Center  (AEDC),  Air  Force  Systems  Command  (AFSC), 

US  Air  Force,  has  been  used  to  extend  previous  experimental  studies  of 
drag  of  spheres  in  low-density  flow.  Using  the  LDH  Tunnel,  drag  data 
have  been  taken  under  hypersonic,  cold-wall  (Tw  «  T0)  support-free  con¬ 
ditions  which  are  described  later. 

These  and  earlier  data  from  other  sources  have  been  examined  in 
regard  to  both  continuum  and  noncontinuum  fluid  flow  behavior.  The  pur¬ 
pose  of  this  report  is  to  present  the  new  data  and  analyses. 

2.0  EXPERIMENTAL  PROCEDURE  AND  DATA 


2.1  LDH  WIND  TUNNEL 

The  L;DH  Tunnel  (Fig.  1)  is  of  the  continuous  type,  with  the  gas 
heated  by  an  electric  arc  or  plasma  generator.  Test  section  conditions 
have  been  determined  from  analysis  of  impact  pressure,  local  mass  flow 
rate,  calorimeter,  and  static  pressure  measurements.  A  description  of 
the  facility  is  contained  in  Ref.  7.  Briefly,  the  tunnel  consists  of  a 
d-c  arc -heater,  stilling  chamber,  conical  nozzle  of  15-deg  half  angle, 
test  chamber  with  instrumentation,  diffuser,  and  pumping  system.  Addi¬ 
tional  information  on  the  nozzle  flow  is  contained  in  Refs.  17  and  18. 

Measurements  of  total  enthalpy  at  the  nozzle  throat  by  calorimetry 
agree  closely  with  total  enthalpy  computed  on  the  basis  of  measured  mass 
flow  rate,  total  pressure,  sonic  throat  area,  and  the  assumption  of  thermo¬ 
dynamic  equilibrium  in  the  fluid  upstream  of  the  throat.  However,  on  the 
basis  that  computed  relaxation  lengths  for  molecular  vibration  downstream 
frcm  the  throat  are  from  10^  to  104  times  local  nozzle  radius,  all  theoret¬ 
ical  evidence  indicates  frozen  flow  from  the  throat  onward.  Thus,  test 
section  flow  characteristics  are  based  on  sudden  freezing  of  the  flow  at  the 
throat. 

Operating  conditions  for  this  experiment  follow: 


Gas 

= 

Nitroge 

n 

Nominal  Nozzle  Station 

- 

1 .  0  in. 

downstream  of  exit 

Dynamic  Pressure,  psfa 

2.  12 

2.  51 

3.  16 

3.31 

Total  Enthalpy,  Btu/'lb 

= 

871 

1110 

1560 

1760 

Total  Temperature,  "R 

= 

3150 

3960 

5400 

6030 

Mach  Number 

= 

10.  8 

10.  6 

10.  7 

10.  6 

Unit  Re^/in. 

= 

466 

426 

344 

301 

Unit  Re2/in. 

= 

47.  0 

42.  6 

37.  6 

34.  7 

Mean  Free  Path,  Aw  .  in. 

= 

0. 0508 

0.0548 

0,  0686 

0.  0776 

using 

- 

1/26  \  y  iLj{ 

Ax,  O 

2 


AEDC-TDR- 62-205 


2.2  SPHERE  DRAG  MEASUREMENT 

2.2.1  Free-Fall  Method 

For  the  smaller  spheres  (0.  0313  to  0.  375-in.  -diam)  a  free-fall 
technique  was  used  to  determine  drag.  The  models  were  injected  into 
the  test  stream  and  their  trajectories  determined  using  a  photographic 
technique  whereby  a  time  exposure  was  made  of  the  moving  model  which 
was  illuminated  by  an  interrupted  light  source.  Spheres  with  diameters 
of  0.  0625  to  0.  375-in.  were  dropped  into  the  test  section  by  an  electro¬ 
magnet  from  a  position  just  above  the  test  section.  However,  the  spheres 
had  a  tendency  to  retain  some  residual  magnetism,  thereby  making  release 
of  the  smaller  spheres  difficult.  In  order  to  facilitate  the  release  of  the 
smaller  spheres,  a  solenoid  was  wound  near  the  end  of  a  short  length  of 
copper  tubing.  With  the  solenoid  energized,  test  spheres  of  0.  0313  to 
0.  156-in.  diam  could  be  suspended  within  the  tubing.  By  deenergizing  the 
solenoid,  the  models  were  free  to  fall  down  the  tubing.  The  end  of  the 
tubing  was  bent  so  that  the  models  could  be  injected  into  the  test  stream 
at  the  appropriate  angle. 

The  light  source  used  to  illuminate  the  spheres  during  their  flight 
through  the  test  section  was  a  1000-watt  projector  lamp.  The  light  inten¬ 
sity  was  modulated  using  a  rotating  disc  containing  holes  drilled  around 
the  periphery.  The  disc  was  rotated  by  a  high-speed  motor  so  that  the 
light  at  the  test  section  could  be  interrupted  at  rates  up  to  1000  per  second. 
By  taking  a  time  exposure  of  the  test  section,  the  trajectory  of  the  sphere 
was  recorded  photographically  at  equal  time  intervals  during  its  flight. 

Since  the  test  section  is  about  1-in.  in  diameter,  a  telephoto  lens  was 
used  so  that  the  photograph  of  the  test  region  could  be  magnified  for  easier 
data  reduction.  An  example  is  shown  in  Fig.  2. 

Surface  temperatures  of  the  spheres  tested  by  the  free-flight  method 
were  approximately  0.  1  total  temperature. 

2.2.2  Axial-Force  Balance 

A  smaller  amount  of  data  was  obtained  earlier  for  spheres  of  0.  265 
to  0.  686-m.  diam  by  means  of  an  external,  sting-type,  water-cooled, 
axial-force  balance  capable  of  measuring  forces  on  the  order  of  0.001  lb. 
Temperatures  of  the  spheres  on  this  balance  were  somewhat  greater  than 
in  the  free-fall  case,  although  Tw  <  T0,  and  minimum  sphere  size  was 
limited  to  approximately  0.  25  in.  A  photograph  of  this  balance  installed 
in  the  LDH  tunnel  is  shown  in  Fig.  3. 
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The  test  spheres  were  steel  ball  bearings  with  a  0.  1-in.  hole  drilled 
tc  admit  the  sting  support.  The  spheres  were  uncooled  except  for  radia¬ 
tion  from  the  surface  and  conduction  through  the  sting  mount.  Wall 
temperature  of  the  spheres  reached  approximately  0.  3  total  temperature. 

It  should  be  noted  that  these  are  essentially  the  same  data  first  pub¬ 
lished  in  Ref.  7,  but  the  points  now  have  been  shifted  slightly  because  of 
the  use  of  a  new  source  (Ref.  19)  for  viscosity  of  nitrogen  at  high  temper¬ 
atures.  Good  agreement  is  evident  in  comparing  free-fall  and  balance 
measurements. 

A  small  quantity  of  new  data  at  higher  Reynolds  numbers  from  the 
50-Inch  Mach  10  Tunnel  (C)  in  VKF  is  included.  These  were  obtained  by 
use  of  the  same  axial-force  balance  used  in  the  LDH  Tunnel. 


2.3  DATA  REDUCTION 

The  data  from  the  free-fall  tests  were  obtained  in  the  form  of  photo¬ 
graphs  (Fig.  2)  which  showed  the  trajectory  of  the  sphere  at  equal  time 
increments.  The  two  larger  light  spots  in  the  example  shown  in  Fig.  2 
were  used  to  orient  the  photograph,  since  the  Polaroid  back  to  the  camera 
could  be  moved  slightly  relative  to  the  tunnel.  These  position  points 
were  obtained  by  a  double  exposure.  One  exposure  was  made  of  the  path 
of  the  sphere  and  the  other  of  the  two  reference  light  sources  which  could 
be  positioned  on  the  nozzle  exit  before  or  after  a  run.  The  position  points 
were  located  approximately  0.  5  in.  above  and  below  the  tunnel  centerline 
and  0.  5  in.  from  the  nozzle  exit. 

The  coordinates  of  each  trajectory  point  on  a  film  were  read  using 
graph  paper  for  a  grid  and  using  the  position  points  for  alignment.  Initially 
it  was  deemed  desirable  to  reduce  the  x,  y  coordinates  obtained  from  the 
photograph  to  an  acceleration  in  the  x-direction  by  using  the  known  accel¬ 
eration  of  gravity  as  a  reference.  Writing  dax/dt*  in  terms  of  the  x,  y  coor¬ 
dinates  with  y  a  function  of  time  gives 

d3x  «**y  A/d_y\  "** 

d'.2  iy  ^di  ^  dy* 

where  d3y/dia  is  the  acceleration  due  to  gravity,  g. 

The  least-squares  method  was  used  to  give  x  *  f(y)  for  polynomials 
of  various  degree,  from  which  dx/dy  and  <I*x/dyJ  were  determined.  The 
term  dy/dt  was  found  by  fitting  a  straight  line  through  the  values  of 
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(y  -  etJ.'2)  plotted  as  a  function  of  t,  i.  e.  ,  if 


then 


Values  of  were  determined  as  a  lunction  of  time  across  the  test 

section.  However,  since  x  =  f(v)  is  not  a  polynomial  for  a  body  in  a 
uniform  force  field,  the  value  obtained  for  iFx/di1  was  not  constant  through¬ 
out  the  test  section.  Therfore,  it  was  concluded  that  a  more  straight¬ 
forward  method  should  be  used. 

Next,  the  difference  in  axial  coordinates  of  adjacent  points,  \x,  was 
plotted  as  shown  in  Fig.  4.  The  slope  of  the  line  through  these  points  is 
proportional  to  the  acceleration  of  the  sphere  in  the  axial  direction.  *  The 
constant  of  proportionality  between  this  slope  and  the  sphere  acceleration 
was  obtained  by  combining  the  known  time  between  successive  points  on 
the  photograph  and  the  known  scale  factor  of  the  photograph.  The  choice 
was  made  to  combine  the  time  and  position  scales  into  one  constant  so  that 
graph  paper  could  be  used  to  read  coordinates  and  the  points  could  be 
plotted  simply  as  a  function  of  point  number. 

Using  Newton’s  second  law  in  the  form 

n  *  m.lVJt*  (2) 

the  total  drag  force  on  a  sphere  was  determined  from  its  acceleration  and 
its  mass.  The  mass  was  determined  by  using  an  analytical  balance  for 
weighing  several  spheres  of  the  same  size  and  then  taking  an  average. 
Values  of  the  drag  coefficient.  On,  were  obtained  for  the  sphere  using  the 
definition 


where  \  is  the  cross-sectional  area  of  the  sphere.  Table  1  summarizes 
the  conditions  and  results  of  these  tests. 

Reduction  of  data  from  the  axial -force  balance  was  accomplished  by 
the  usual,  well-known  procedure  involving  use  of  a  calibration  curve 
obtained  by  static  calibration  of  the  balance  using  known  drag  loads. 


*It  can  be  shown  that  this  finite  difference  method  for  obtaining  the 
acceleration  of  a  model  in  a  uniform  force  field  gives  a  result  equal  to 
that  from  the  second  derivative  of  the  position-time  curve. 
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3.0  ANALYSIS 


3.1  THE  SPHERE  IN  LOW-DENSITY,  CONTINUUM  FLOW 

It  is  a  generally  accepted  approximation  that  the  drag  coefficient  of 
a  sphere  in  supersonic,  low-density  but  continuum  flow  may  be  expressed 
as  the  sum  of  its  drag  coefficient  at  essentially  infinite  Reynolds  number, 
C'nj ,  plus  the  contribution  from  skin  friction,  plus  a  third  term  represent¬ 
ing  the  combined  influences  of  vorticity,  slip,  temperature  jump,  and 
non- negligible  ratios  of  boundary- layer  thickness  to  sphere  radius,  cf 
Refs.  11,  20-22. 


Consideration  of  only  the  dominant  inviscid  and  first-order  friction 
terms  led  Chahine  (Ref.  13)  to  the  result  (in  the  present  nomenclature). 


Cn  -  Cj){ 


9  C/8 

Rea  P~/P2 


(4) 


where  C"  is  the  shear  function  given  by  Cohen  and  Reshotko  in  Ref.  23. 
The  9/8  is  an  empirically  derived  constant  apparently  based  on  data  taken 
under  conditions  where  M*,  =4,  Tw  2:  t,,  and  Ta  -  T0.  Chahine  treats  Fw 
as  a  constant,  taking  the  stagnation  point  value. 


Rott  and  Wittenbury  in  Ref.  11  considered  all  the  factors  mentioned 
above  and  arrived  at  the  expression  (again  using  present  nomenclature): 

r/iptV„  Bind)  ~  Ct  V(VR)  U3/UJ  +  C2  ( *2/R  )  (a/U^, )  +  ■  •  •  (5) 

where 


Ct  -  f.  lPo*/p2.  /*(T),  Pr,  . . 1 

C2  —  f2  (vorticity,  displacement,  curvature,  slip,  •  •  ■  •  ) 


The  form  of  Eq.  (5)  may  be  made  more  nearly  similar  to  that  of  Eq.  (4) 
by  using  the  following  relations: 

■°oon™  53  P2U2 

H/Aj  =  l/Krij  *  const.  Rej/Mj 

=  const,  ipfx/p R/(Aj  Re2) 


Substituting  the  above  relations  into  Eq.  (5)  yields 


. ■  -  Z  sin  ■>  : 


1 
•  I 

J 
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If  K,  is  regarded  as  a  function  of  ,  among  other  things,  then  one  is 

led  to  write 


Cn  -  Cn,  4  - .  (6) 

v  Re,  Re, 

This  is  the  form  indicated  by  Aroesty  (Ref.  16)  in  discussing  the  subject, 
and  it  is  interesting  to  note  his  comment  that  neither  the  sign  nor  mag¬ 
nitude  of  K,  has  been  established  previously.  The  "approximately  equals” 
symbol  is  inserted  by  the  present  authors,  partly  out  of  consideration  of 
the  finite  range  of  validity  of  the  equation  and  partly  because  the  coeffi¬ 
cients  K,  and  K,  are  later  determined  empirically. 

A  presentation  of  both  the  new  and  the  more  recent  of  the  previously 
published  measurements  of  drag  of  spheres  in  supersonic  and  hypersonic 
flows  is  given  in  Figs.  5-7  where  Re,  is  used  as  the  parameter.  *  Curves 
of  the  form  described  by  Eq.  (6)  have  been  fitted  to  these  data,  with  Cn. 
taken  from  Ref.  24  for  corresponding  Mach  numbers.  The  values  of  K, 
and  K,  determined  by  this  process  are  presented  in  Figs.  8  and  9  as 
functions  of  the  enthalpy  parameter, 

Hw  -  hw/ho  -  1  (7) 

because  the  heat-transfer  situation  would  be  expected  to  be  important. 
However,  it  is  apparent  that  Mach  number  or  some  function  thereof  exerts 
an  influence  on  the  coefficients. 

Although  there  is  the  possibility  of  correlating  K,  and  possibly  K,  with 
known  parameters,  it  must  be  remembered  that  base  drag  could  be  impor¬ 
tant  at  the  lower  Mach  numbers  (Refs.  25-26).  In  view  of  this,  plots  of  the 
coefficients  as  functions  of  Mach  number  or  density  ratio  may  not  be  unique 
for  Mm  <  6.  Therefore,  no  adc  itional  correlations  are  presented.  There  is 
evidence  in  Figs.  8  and  9  that  <,  and  K,  tend  toward  dependence  on  Hw 
alone  at  large  Mach  numbers  in  the  essentially  perfect  gases  represented. 
One  would  expect  this,  excepting  effects  such  as  chemical  relaxation  and 
ablation  which  were  not  present  in  the  tests  being  discussed. 

It  is  thought  that  the  more  conclusive  portions  of  Figs.  8  and  9  are 
those  pertaining  to  the  cold- wall,  hypersonic  case  of  >  6  and  Hw  ■*  -1, 
because  the  influence  of  base  pressure  and  other  Mach  number  dependent 
factors  must  be  minimized  there.  Also,  variations  about  the  nominal 


*The  conventional  thin  shock  wave  and  inviscid  shock  layer  are 
assumed  in  computing  Re,,  and  the  sphere  diameter  is  used  as  the  charac¬ 
teristic  length. 
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Mach  number  in  the  experiments  would  have  less  effect  on  the  data, 
thereby  decreasing  scatter  and  aiding  in  determining  the  coefficients. 

The  situation  is  less  clear  in  regard  to  Ka  because  it  assumes  im¬ 
portance  only  at  the  lowest  Reynolds  numbers  and  some  older  data  do 
not  penetrate  sufficiently  into  that  regime.  The  values  of  K2  deduced 
from  the  data  are  consistently  negative,  which  is  rather  interesting 
because,  to  the  authors'  knowledge,  this  possibility  has  not  been  sug¬ 
gested  previously. 

Van  Dyke  (Ref.  27)  has  calculated  an  example  pertaining  to  the  stag¬ 
nation  region  of  a  sphere  with  entropy  gradient,  longitudinal  curvature, 
transverse  curvature,  slip,  and  temperature  jump  considered.  He 
assumed  that  M«,  =  Pr  «  0.7,  p  ~  T,  Hw  =  -0.8,  and  y  «  7/5,  and  neg¬ 
lected  the  displacement  boundary  laye.^.  One  result,  which  may  be 
expressed  in  terms  of  and  K* ,  is,  in  the  present  nomenclature, 

Kj/K,  {at  stag.  pt. )  =  .  0.661 

On  the  other  hand,  we  find  experimentally,  for  Hw  =  -0.8  and  Mw  -  1« , 

Kj/Kj.  (average)  —  —  0.57 

where  the  negative  sign  is  due  tc  (average).  (The  sign  of  K,  causes 
to  be  larger  than  it  would  be  if  K2  were  neglected  as  in  Eq.  (4). 

It  should  be  noted  that  the  terms  in  the  algebraic  sum  giving  K2/Kt  from 
Van  Dyke's  analysis  are  both  positive  and  negative  in  sign.  Therefore,  it 
is  not  inconceivable  that  K,/K,  (average)  would  be  less  than  the  stagnation 
region  value  or  that  it  would  be  negative.  It  may  be  relevant  to  add  Van 
Dyke's  remark  that  the  effect  of  entropy  gradient  calculated  by  Van  Dyke  is 
twice  as  large  as  that  found  by  Probstein  (Ref.  28),  but  less  than  one-half 
that  predicted  by  Ferri,  Zakkay,  and  Ting  (Ref.  29).  Also,  the  Navier- 
Stokes  equations  on  which  Van  Dyke's  analysis  is  based  may  not  be  valid 
below  Re2  =  0(100),  but  the  experimentally  determined  coefficients  were 
strongly  influenced  by  data  corresponding  to  Re*  «  100. 

In  view  of  the  uncertainty  that  would  exist  in  theoretical  estimates  of 
the  coefficient  Ka  (average),  it  appears  that  the  experimental  values  pre¬ 
sented  here  are  reasonable.  No  other  experimentally  derived  values  are 
known  to  the  writers. 

It  is  important  to  keep  in  mind  that  Eq.  (6)  must  fail  at  sufficiently  low 
Reynolds  numbers.  This  failure  may  be  postponed  by  inserting  additional 
terms,  but  that  is  not  possible  in  the  absence  of  reliable  drag  data  for  very 
nearly  free- molecular  flow'.  Finding  the  subsequent  terms  does  not  seem 
highly  important  in  the  cold-wall,  hypersonic  case  since  the  two  terms 
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identified  here  are  sufficient  to  the  limit  of  continuum  flow  where  it  is 
then  necessary  to  seek  further  extensions  by  means  of  noncontinuum  flow 
analysis. 


3.2  THE  SPHERE  IN  NONCONTINUUM  FLOW 


3.2.1  Some  Comments  on  Mean  Free  Path 


Often  in  the  study  of  elementary  kinetic  theory  of  gases  the  gas  mole¬ 
cules  are  approximated  by  small,  perfectly  elastic  spheres  (billiard-ball 
model).  This  simplified  model  naturally  leads  to  the  idea  of  a  molecular 
collision  and  suggests  the  concept  of  a  free  path  between  collisions.  The 
mean  free  path,  A,  of  these  billiard-ball  molecules  can  be  expressed  as 

K  -  no)  (8) 

where  n  is  the  number  density  of  molecules  and  o  is  the  diameter  of  the 
molecule. 


This  equation  is  of  little  use  for  calculating  mean  free  path  since  o 
cannot  be  determined  directly.  However,  the  equation  shows  that  the 
mean  free  path  for  an  elastic  sphere  molecule  is  a  function  of  gas  density 
only  and  is  inversely  proportional  to  it. 


The  idea  of  the  billiard-ball  molecule  can  be  extended,  and  expres¬ 
sions  for  the  transport  properties  of  viscosity  and  thermal  conductivity 
can  be  obtained.  These  expressions  are,  respectively, 


where  v  is  the  mean  molecular  speed,  P  the  mass  per  unit  volume,  and 
Cv  the  specific  heat  at  constant  volume.  The  most  simplified  calculation 
gives  the  value  of  1/3  for  both  constants  of  proportionality. 


More  comprehensive  analyses  have  yielded  values  for  the  constant 
of  proportionality  such  as  0.350,  0,  310,  0.  491,  and  0.  499  in  the  formula 
for  viscosity.  Values  for  the  constant  of  proportionality  in  the  formula 
for  thermal  conductivity  vary  from  1.  10  to  2.  52  times  greater  than  those 
for  viscosity. 


The  many  different  constants  obtainable  make  the  numerical  deter¬ 
mination  of  mean  free  path  very  questionable.  The  only  encouragement 
comes  from  the  fact  that  all  analyses  give  the  same  general  form  for  the 
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equations.  Other  models  for  the  molecule  besides  the  hard  sphere  also 
give  equations  of  this  form.  Any  assumed  force  field  between  molecules 
except  the  billiard-ball  type  does  away  with  the  idea  of  a  rigorous  free 
path.  Each  molecule  is  in  contact  with  all  other  molecules,  to  some 
extent,  at  all  times. 

The  equations  mentioned  above  can  be  applied  to  obtain  an  effective 
mean  free  path  for  each  property.  Thus,  there  is  obtained  an  effective 
mean  free  path  for  viscosity  and  one  for  thermal  conductivity.  These 
two  mean  free  paths  are  different  and  vary  differently  with  temperature. 
These  differences  can  be  explained  on  the  grounds  that  distinctions  in 
different  molecular  properties  are  important  in  defining  different  effec¬ 
tive  mean  free  paths.  Viscosity  is  associated  with  the  transfer  of 
momentum  across  a  plane,  and  the  exchange  of  momentum  during  a  col¬ 
lision  is  important,  whereas  thermal  conductivity  arises  from  the 
exchange  of  energy  during  a  molecular  encounter.  Each  of  the  mentioned 
phenomena  is  associated  with  a  central  type  collision.  (A  collision  in 
which  the  molecules  are  only  slightly  deflected  does  not  exchange  an 
appreciable  amount  of  momentum  or  energy.  )  Other  phenomena,  such  as 
scattering,  would  place  just  as  much  importance  upon  the  glancing  col¬ 
lisions.  Molecular  scattering  experiments  give  mean  free  paths  that  are 
less  than  half  those  from  viscosity  measurements. 

These  effective  mean  free  paths  can  be  used  with  Eq.  (8)  to  define 
molecule  diameter  or  collision  cross  sections.  It  is  found  that  the 
collision  cross  section  decreases  with  rising  temperature.  These  results 
are  expected  for  molecules  with  a  force  field  different  from  the  billiard- 
ball  model.  The  effect  of  an  increase  in  temperature  is  to  increase  the 
molecular  velocities.  A  molecule  will  therefore  reside  a  shorter  time 
in  the  vicinity  of  another  molecule  and  experience  a  correspondingly 
smaller  deflection  which  is  equivalent  to  a  smaller  effective  collision  cross 
section. 

In  order  to  bypass  these  problems,  the  viscosity  equation  may  be  used 
to  determine  the  form  for  the  mean  free  path  as  a  function  of  temperature, 
but  the  constant  will  be  unspecified.  In  other  words,  instead  of  using  mean 
free  path  as  a  parameter,  iWy/ipa)  may  be  used.  The  mean  molecular 
velocity,  v,  has  been  replaced  by  the  speed  of  sound,  «.  This  may  be  done 
because  a  is  directly  proportional  to  v  for  constant  ratio  of  specific  heats, 
y,  and  is  given  by  the  following  equation: 


The  choice  to  replace  v  by  a  may  be  justified  because  the  speed  of  sound 
appears  more  frequently  in  the  study  of  high-speed  gas  flows,  and  tabu¬ 
lated  values  of  the  speed  of  sound  are  readily  available  for  real  gases. 
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The  dimension  of  nyy/ipa)  is  the  same  as  that  of  mean  free  path, 

1.  e.  length.  This  parameter  can  be  ratioed  to  a  body  dimension  to 
obtain  a  dimensionless  ratio.  Using  the  sphere  diameter,  d,  for  this 
characteristic  dimension,  the  dimensionless  parameter  becomes 
pVy/fdpa).  Multiplying  both  numerator  and  denominator  of  this  expres¬ 
sion  by  the  gas  velocity  gives 


Kn  -  v7  y  M/Re  (9) 

It  has  been  shown  that  the  mean  free  path  of  real  gas  molecules  is 
not  uniquely  determined.  However,  in  subsequent  theoretical  develop¬ 
ments  it  is  assumed  that  some  effective  mean  free  path  exists  and  that 
the  gas  molecules  can  be  treated  as  elastic  spheres.  In  order  to  com¬ 
pare  a  theoretically  derived  equation  with  experimental  data,  a  free  con¬ 
stant  is  later  introduced  in  the  relationship  between  mean  free  path  and 
the  quantity. 


1 1  \/  y  /(pa) 


(10) 


3.2.2  Interaction  of  Gas  Molecules  and  a  Solid  Surface 


One  of  the  most  difficult  problems  arising  in  a  kinetic  theory  dis¬ 
cussion  of  gas  flo  is  that  of  the  interactions  between  the  gas  molecules 
and  a  solid  boundary.  The  problem  of  this  gas-solid  interaction  is 
associated  with  the  reflection  or  reemission  of  the  molecules  from  the 
surface,  i.  e.  ,  how  the  momentum  and  energy  of  the  molecules  and  their 
different  degrees  of  freedom  are  accommodated  to  wall  conditions.  For 
polyatomic  molecules  there  is  an  accommodation  for  each  mode  of  inter¬ 
action.  Thus,  there  is  an  accommodation  coefficient  of  mass,  transla¬ 
tional  energy,  vibrational  energy,  tangential  momentum,  normal  momentum, 
etc.  The  reader  interested  in  a  detailed  discussion  of  the  interaction  phe¬ 
nomenon  may  find  several  reviews  of  the  subject  (cf  Ref.  30). 

Experimental  results  indicate  that  for  usual  engineering  surfaces  the 
actual  reflection  process  is  near  a  diffuse  reflection.  The  equation 
describing  a  diffuse  molecular  reflection  with  complete  accommodation  to 
wall  conditions  will  be  developed  below. 

The  number  of  molecules  incident  upon  an  elemental  area,  dA,  per  unit 
time  is,  from  Ref.  31  (See  nomenclature  and  Fig.  10.), 

n,  /  k  1  e  cos  ^ 1  +  v'ff  cos  0)  •  1  +  erf  (S  cos  0)  1  1  d  A  (11) 

For  the  development  of  diffuse  reflection  from  a  surface,  it  is 
assumed  that  the  reflection  from  each  incremental  surface  area  is  the 
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same  as  the  effusive  flow  of  a  hypothetical  gas  at  wall  temperature 
through  an  opening  replacing  the  incremental  area.  The  density  of  this 
hypothetical  gas  can  be  determined  from  the  conservation  of  molecules. 


Kennard  (Ref.  32)  gives  the  number  of  molecules  crossing  any  incre¬ 
mental  area  within  a  gas  per  unit  time  with  speed  between  v  and  v  j-  dv 
within  the  solid  angle,  dej,  as 


a  w 


(  2  rr  kTw) 


9  dv  d  <i>  dA 


(12) 


In  the  free- molecule  regime  the  flow  of  molecules  through  an  opening, 
dA,  is  the  same  as  that  given  in  Eq.  (12).  Integrating  over  all  speeds,  the 
total  number  of  molecules  passing  (reemitted  from)  dA  within  the  solid 
angle,  dw.  is 


(  2  tt  kT ) 


cos  <f}  do> 


dAX 


(^f  ( 


cos  ocj  d4\(13) 


Because  it  is  assumed  that  the  flow  is  reflected  symmetrically  about 
the  normal  to  the  surface,  do >  may  be  written  as  (See  Fig.  10.  ) 

da>  =  2^sin  d  4>  (14) 

The  total  number  of  molecules  within  the  angle  <i>  to  <p  +  d<£  is, 
substituting  Eq.  (14)  in  Eq.  (13), 

„w  ±±±  JA  (15) 

'  m  '  yj  n 


Integrating  Eq.  (15)  over  all  directions  in  which  molecules  can  be 
reflected  from  the  surface  (o  ^  gives 


(  JJlLlV  2  d  A 

\  nm  f  i 


(16) 


the  total  number  of  molecules  reflected  from  the  surface.  Equating  this 
expression  to  the  number  of  incident  molecules  (Eq.  (11))  gives  a  solution 
for  nw: 


nw 


-!jf-  je  (  5  co 8  ^  -i ■  \/  n  (S  cos  9)  ^  l  +  erf  (S  cos  0 )  !  |  (17) 
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The  mean  speed  of  the  reemitted  molecules  is  defined  as 

/  *  s/ a  yc  “  j  tT  ~~ 

“w  )  cosodcydA  v4  e  w  dv 

vw  - - - - 1 

l  »  \2  .  ,  ,  f  00  ,  "UT»  j 

nw  [ =— —  cos  ?!  do  dA  v  e  dv 

w  \2  TTtTy,/  *  Jo 


or 


kT 

m 


(18) 


3.2.3  Derivation  of  Mean  Free  Path  of  Reemitted  Molecules 

In  any  discussion  of  free-molecular  flow  or  near-free-molecular 
flow,  the  dimensionless  ratio  of  some  mean  path  of  gas  molecules  to 
some  characteristic  dimension  usually  is  introduced  as  a  parameter. 

This  dimensionless  ratio  is,  of  course,  the  Knudsen  number.  One  prob¬ 
lem  that  arises  in  connection  with  the  Knudsen  number  is  that  of  deter¬ 
mining  which  mean  free  path  is  the  most  significant  in  the  problem  at 
hand.  Another  important  related  question  is  that  of  the  numerical  deter¬ 
mination  of  this  mean  free  path.  This  latter  problem  was  discussed  in 
Section  3.2.1. 

The  Knudsen  number  may  be  based  upon  a  "free-stream"  mean  free 
path  which  is  independent  of  free-stream  velocity  and  body  conditions. 
However,  if  interference  between  free-stream  and  reemitted  molecules 
exists,  it  is  reasonable  to  seek  some  mean  free  path  more  descriptive  of 
the  flow.  Since  the  molecular  mean  free  path  is  always  finite,  the  mole¬ 
cules  reflected  from  the  body  will  always  disturb  the  free  stream  to  some 
extent.  A  criterion  for  closeness  to  truly  free-molecule  flow  is  the 
average  distance  a  reflected  molecule  travels  after  leaving  the  body  sur¬ 
face  before  experiencing  a  collision  with  a  free-stream  molecule.  This 
mean  free  path  of  reemitted  molecules,  A.w,  is  determined  in  the  remain¬ 
der  of  this  section. 

The  velocity,  relative  to  the  free  stream,  of  a  molecule  leaving  the 
body  is  the  vector  sum,  v  +  UM.  The  mean  free  path  of  a  reflected  mole¬ 
cule  is  different  from  A^  since  the  velocity  is  different  from  the  mean 
thermal  velocity  in  the  free  stream,  ¥M.  Jeans  (Ref.  33)  gives  the  mean 
free  path  of  molecules  with  speed,  v,  in  a  static  gas.  This  mean  free 
path  for  a  molecule  with  speed,  «,  is 

A  (v)  =  Am  Gtv/v,.)  (19) 

A  graphical  representation  of  this  is  given  in  Fig.  11.  (The  equation  of 
G(v/vw)  is  given  later  as  Eq.  (43).  ) 
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Therefore,  the  distance  traveled  by  the  reflected  molecules  is 

A  (  ,  »  -  G  ^ 

Notice  that  this  is  not  the  distance  from  the  body,  because  the  body  has 
been  moving  with  velocity  since  the  molecule  left  the  surface.  The 
time  required  to  traverse  this  free  path  is 

i^UMl 

Converting  to  a  coordinate  system  fixed  relative  to  the  body  and  using  the 
velocity  of  the  reflected  molecule  and  the  mean  time  before  a  collision  as 
developed  above,  the  mean  free  path,  relative  to  the  body,  of  a  molecule 
with  velocity  v  is 


A.  G 

l  v  +  O 


(20) 


If  it  is  assumed  that  the  directional  distribution  of  mean  free  paths 
is  similar  for  spheres  under  varying  conditions,  then  a  representative 
free  path  that  is  easy  to  calculate  can  be  chosen.  The  mean  free  path 
easiest  to  determine  and  the  most  important  from  the  standpoint  of  alter¬ 
ing  the  free  stream  is  the  free  path  of  molecules  directed  upstream  into 
the  oncoming  flow.  Expressed  in  terms  of  the  average  speed  of  the  mole¬ 
cules  leaving  the  surface  this  reemitted  mean  free  path*,  A.w ,  becomes 


(21) 


The  mean  velocity  of  molecules  leaving  the  sphere  is,  from  Eq.  (18), 


V 


--  R  T* 


*Probstein  (Ref.  34)  defines  a  mean  free  path  of  reemitted  molecules 
relative  to  incident  molecules  for  a  hypersonic  Mach  number  and  a  highly 
cooled  body.  He  concludes  that  this  mean  free  path  is  the  mean  free  path 
most  significant  in  the  free-molecule  flow  and  near-free-molecular  regimes. 
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The  free-stream  velocity  may  be  written  as 

(22) 

and  the  free-stream,  mean  thermal  velocity  is 


8RTm 


(23) 


Thus,  substituting  Eqs.  (18),  (22),  and  (23)  into  Eq.  (21)  gives  the  mean 
free  path  of  molecules  reflected  from  the  body  with  velocity  vw  and 
directed  upstream: 


/iiL  M  iyn-) 
v 6  v  T~  v  8  / 


*_L  (1 

9ff  \  T 


(24) 


m 


Dividing  both  sides  of  Eq.  (24)  by  a  reference  length,  one  finds 


Knw  =  Kn„ 


(W: 


j-Ll.) 

\f  8  / 


oy  /  »  i 
9  n  \  Tw  > 


(25) 


Equation  (25)  is  based  upon  the  assumption  of  complete  accommoda¬ 
tion  of  the  reemitted  moleculeB  to  the  wall  condition.  For  specular 
reflection  from  the  surface,  vw  is  equal  to  the  incident  velocity,  ,  and 
Eq.  (21)  becomes 

A.  -  A«-  G  (26) 

Independent  of  the  type  of  molecule  reflection,  the  reflected  mean  free 
path  is  proportional  to  the  free-stream  mean  free  path  for  constant  free- 
stream  conditions  and  wall  temperature,  i.  e.  , 

Kn...  ^  Kn  for  M_  **  const,  and 
T^/Tw  =  const. 

Combining  this  relationship  with  Eq.  (9)  gives 

Knw  -  VF  <27> 

The  mean  free  path  of  molecules  emitted  from  the  body  and  directed 
upstream  will,  in  general,  not  equal  the  mean  free  path  of  molecules  with 
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the  mean  emitted  velocity,  vw  .  The  rigorous  definition  of  mean  free 
path  is 

G  S) 

fo  (y  /  c  )  v*atJdv  (28) 

/  n  (v)  v*  d  <o  dv 


where  n(v)  is  the  number  density  of  reemitted  molecules  in  velocity- 
position  coordinates  and  is  given  by  Eq.  (12).  The  term  dtuv3  dv  is  the 
elemental  volume  in  velocity  space. 

For  hypersonic  flows,  »  l.  Therefore,  from  Fig.  11, 

V 


If  the  body  is  highly  cooled,  v  «  over  most  of  the  range  of  v  and  the 

term  (v  +  U^)  can  be  replaced  by  .  Evaluating  the  mean  free  path  of 

reflected  molecules  for  highly  cooled  bodies  in  hypersonic  flow  gives, 
from  Eq.  (28), 


.  A  in  ils  (29) 

00  Lw  V  8  m 

This  result  is  the  same  as  that  obtained  from  Eq.  (24)  after  assuming  the 
conditions  of  hypersonic  flow  and  a  highly  cooled  body.  For  conditions 
other  than  the  highly  cooled  body  in  a  hypersonic  stream  it  will  be  assumed 
that  the  mean  free  path  of  reemitted  molecules  as  given  by  Eq.  (24)  is 
applicable. 


3.2.4  Aerodynamic  Drag  in  Near-Free^Molocular  Flow 

In  the  near-free-molecular  flow  regime  the  number  of  free- stream 
molecules  striking  the  surface  is  diminished  from  the  free-molecular 
value  because  some  have  been  deflected  by  molecules  reemitted  from  the 
body.  In  order  to  facilitate  a  discussion  of  molecular  collisions,  velocity 
distribution  of  molecules,  and  other  properties  of  molecules,  it  is  desir¬ 
able  to  divide  the  molecules  into  various  classes  according  to  their  past 
history.  The  undisturbed  medium  through  which  the  body  is  moving  is 
identified  by  the  subscript  «  .  Molecules  that  have  been  reflected  from 
the  wall  and  have  not  experienced  a  collision  with  molecules  of  a  different 
class  are  given  the  subscript  w  .  Molecules  that  actually  strike  the  sur¬ 
face  are  known  as  incident  molecules,  i  ■  In  the  near-free-molecule 
region  it  will  be  assumed  that  the  incident  molecules  are  composed  of 
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molecules  from  the  free  stream  that  have  not  undergone  a  collision  with 
molecules  emitted  from  the  wall  prior  to  impacting  upon  the  surface. 
Free- stream  molecules  and  molecules  from  the  wall  that  collide  with 
each  other  are  shifted  into  the  class  of  scattered  molecules,  s  . 


In  the  near-free-molecule  flow  analysis  that  follows,  it  is  assumed 
that  a  free- stream  molecule  is  put  in  the  scattered  class  only  by  a  col¬ 
lision  with  a  molecule  from  the  wall,  and  that  scattered  molecules  do  not 
contribute  to  the  drag  of  the  sphere.  In  other  words,  the  shielding  effect 
of  the  molecules  emitted  from  the  sphere  is  the  only  effect  considered. 

Two  analytical  approaches  have  been  attempted  for  the  near-free- 
molecule  flow  regime.  One  of  these  i6  an  intuitive  approach  wherein  the 
effects  of  collisions  between  molecules  emitted  from  the  body  and  free- 
stream  molecules  are  considered.  This  is  the  so-called  first-collision 
method.  The  other  approach  is  an  attempt  to  approximate  the  solution  of 
the  Boltzmann  equation  for  a  body  in  near-free -molecule  flow.  Among  tne 
noteworthy  investigations  are  those  of  Baker  and  Charwat  (Ref.  9)  who 
followed  the  first-collision  method  and  that  of  Willis  (Ref.  10)  who  obtained 
an  iterative  solution  of  the  Boltzmann  equation.  A  semiempirical  method 
for  the  transitional  region  has  been  proposed  by  Rott  and  Whittenbury 
(Ref.  11).  The  results  of  these  studies  will  be  discussed  briefly  below. 

Baker  and  Charwat  (Ref.  9)  consider  the  first  collision  of  molecules 
reflected  from  a  sphere  with  free- stream  molecules.  The  equation  they 
obtained  for  Co  is 


Cn  -  2 


1  +  h.  (0.444)  -  B  ( 0.23 )  -  B(0.65) 

Vo  vo 


(30) 


where  terms  of  order  Ba  have  been  neglected.  Retaining  these  terms, 
their  analysis  gives 


(31) 


(,  ♦  i 

(0.444)  -  B 

r  -  j 

1  0.386  (—-)  +  0.637 

+  0.23 

1 

|  V  v0/ 

V. 

+  B*  1 

0.566 

+  0.181  -  0.007 

1 

V0 

J 

where  B  and  Ve/V0  are,  as  calculated  by  the  present  authors, 
B  =  1/Knw 

and 


v _  _Ve_  =  f~2n_  /1\  iTw 

“  V.  -  4  \  y  V  Tw 


(32) 
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This  expression  for  Ve/V0  is  obtained  if  Vc  is  assumed  to  be  the  average 
velocity  of  molecules  passing  a  plane  in  a  gas  at  a  temperature  of  Tw. 
This  value  of  Ve/V0  gives  the  correct  free- molecular  value  of  Cn  for  the 
case  of  B  -*  0  .  The  value  obtained  in  Baker  and  Charwat’s  paper  is 


v0 


*  (7)  (x>  (ft)' 


which  apparently  is  based  upon  a  value  for  Ve  of  v"  3  RT w  or  the  root  mean 
square  of  the  velocity  of  gas  molecules  contained  in  a  volume  of  gas  at 
a  temperature  equal  to  Tw .  It  is  necessary  to  use  Eq.  (32)  in  order  to 
obtain  the  correct  free-molecule  limit. 


Willis  (Ref.  10)  obtained  a  near-free-molecular  approximation  to 
the  Boltzmann  equation  and  applied  it  to  obtain  the  drag  coefficient  of  a 
sphere.  He  considered  the  solution  for  both  hard  sphere  molecules  and 
Maxwellian  molecules  and  found,  for  Maxwell  molecules. 


Cd  ■  2 


0.59 

Sh 


°  (°-17  *  -it) 


and,  for  hard  sphere  molecules. 


Cn  *  2 


0.59 


a  (0.3?  + 


(33) 


(34) 


In  terms  of  this  paper. 


and 


Sb  «*  s/y/2 


y/y 

4  Kaw  yf  Tw 


(35) 


This  value  of  a  is  appropriate  only  for  the  hard  sphere  molecular  model. 
However,  for  comparison,  this  definition  of  a  will  be  used  in  the  drag 
equation  for  the  Maxwell  molecules. 


For  hypersonic  flow  Eq.  (25)  may  be  simplified  to 


Knw 


nr 


k) 


<36) 


Equation  (36)  may  be  used  with  Eq.  (35)  to  express  a  in  terms  of  Knw. 
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Rott  and  Whittenbury  {Ref.  11)  assume  that  the  flow  in  the  transi¬ 
tional  region  can  be  represented  by  a  two-component  fluid,  one  com¬ 
ponent  is  treated  as  a  molecular  beam,  whereas  the  other  component 
is  treated  as  a  continuum  fluid.  Their  method  is  semiempirical  in  the 
sense  that  the  two  extremes,  i.  e.  ,  the  free  molecular  and  the  inviscid 
value  of  Co,  are  obtained  from  experiment.  They  obtain  an  expression 
for  Cg  of  the  form 


Cd  =  CDi  +  e  C>/Kn»  (CDfin  -  Cdi)  < 37) 

where  C,  is  a  "free"  constant  to  make  the  equation  best  fit  the  experi¬ 
mental  data. 

In  a  discussion  of  kinetic  theory  it  is  assumed  that  the  dynamics  of 
a  molecule  are  independent  of  its  past  history.  To  be  more  specific,  it 
is  assumed  that  the  probability  of  a  molecule  undergoing  a  collision  within 
the  next  instant  of  time  is  not  a  function  of  the  time  since  the  last  collision. 
Let  a  (v)  d  ft  be  the  probability  that  a  molecule  moving  with  velocity  v  in 
the  direction  of  ft  will  collide  with  other  molecules  within  the  distance  df  , 

Instead  of  looking  at  only  one  molecule  and  giving  odds  on  its  being 
in  a  collision,  we  can  consider  a  large  number  of  molecules  and  deter¬ 
mine  the  distribution  of  collisions  with  other  molecules  as  a  function  of 
length.  Consider  a  large  number  of  molecules,  N(0),  leaving  a  point, 
ft  =  0,  with  velocity,  v,  in  the  direction  of  a  coordinate  axis,  ft.  The  N(0) 
molecules  are  emitted  at  intervals  of  time  large  enough  to  prevent  inter¬ 
ference  with  each  other.  N(i)  is  the  number  of  the  emitted  group  of  mole¬ 
cules  at  a  distance  ft  from  the  initial  point  that  have  not  experienced  a 
collision. 

The  number  of  molecules  that  undergo  a  collision  between  ft  and 
4+  &  ft  can  be  expressed  as  the  product  of  the  number  of  molecules  at  it 
and  the  probability  that  a  single  molecule  will  have  a  collision  in  Aft. 

This  may  be  written  as 

N  (I)  a  (v)  d  jC  (38) 

The  number  of  colliding  molecules  can  also  be  regarded  as  the  negative 
of  the  rate  of  change  of  along  the  direction  of  motion  times  the  incre¬ 
mental  length,  i.  e.  , 


-  r  dN  1  dJ0  (39) 

Equating  quantities  (38)  and  (39)  yields  the  relation 

-  [dN(jg)/d|]di  =  N  U)  a  (v)d  £  (40) 
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Solving  this  differential  equation  and  applying  the  boundary  condition 
that  N(i)  =  N  (0)  at  J?  =  0  gives 

N  U)  =  N  CO)  e  (vU  <41> 

The  mean  free  path  of  molecules  with  speed,  vt  is 


which  reduces  to 


A  (v)  «  1 la  (v) 


(42) 


Jeans  (Ref.  33)  gives  the  mean  free  path  of  gas  molecules  as  a 
function  of  their  velocity  relative  to  a  static  gas.  His  result  is  presented 
as  a  curve  in  Fig.  11.  The  equation  of  this  curve  is 


6fc)' 


A  (v 
A 


8  V  2 


M. 


a 


•*[££)]} 


(43) 


It  is  assumed  that  the  mean  free  path  of  molecules  reflected  from  a 
wall  is  equal  to  the  mean  free  path  of  a  molecule  with  velocity  vw,  where 
»w  is  the  mean  velocity  of  molecules  emitted  from  the  wall.  This  mean 
free  path  will  be  designated  Aw .  The  number  of  molecules  that  experience 
a  collision  between  distances  of  i  and  S.  +  di  from  the  wall  is  therefore 

N  (0)  e“  iMw  (44) 

This  expression  will  be  used  to  obtain  the  total  number  of  collisions 
between  reemitted  and  free-stream  molecules  at  a  point  in  space. 

If  the  mean  free  path  of  gas  molecules  reflected  from  the  surface  is 
very  large  compared  to  the  sphere  diameter,  but  still  finite,  a  molecule 
leaving  the  model  can  be  considered  to  originate  from  a  point  source.  A 
simplified  analysis  will  follow  in  which  it  is  assumed  that  the  sphere  can 
be  considered  a  point  source.  A  more  comprehensive  analysis  is  under¬ 
taken  in  the  next  section  wherein  the  finite  size  of  the  sphere  is  taken  into 
account. 

By  determining  the  relative  number  of  collisions  per  unit  time  be¬ 
tween  free-stream  molecules  and  molecules  from  the  wall,  the  number- 
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of  free- stream  molecules  that  become  incident  molecules  can  be  deter¬ 
mined.  Since  it  is  assumed  that  incident  molecules  are  composed  of 
only  free- stream  molecules  that  have  not  experienced  a  collision,  the 
drag  coefficient  of  the  sphere  is  simply 

cd  -  CDfm  (45) 

in  the  near-free-molecule  regime. 

Consider  the  molecules  leaving  a  surface,  dA,  within  a  solid  angle, 
d  to.  These  molecules  will  remain  within  this  solid  angle  until  they  under¬ 
go  a  collision  with  another  molecule.  The  distribution  of  collisions  of 
these  N  f  0)  molecules  is  given  by  Eq.  (43).  The  number  of  molecules 
within  d<u  must  be  determined  by  an  integration  over  the  sphere  surface. 

As  developed  earlier,  the  number  of  molecules  from  each  element 
of  surface  area,  dA ,  within  dw  is,  from  Eq.  (13), 

i/a 

/  2  kT„\  cos  (f>  dot  dA 
'  ra  '  2  U)i/2 

where  nw  is  given  by  Eq.  (17).  Confining  attention  to  the  reflected  mole¬ 
cules  which  are  directed  upstream  into  the  incident  flow,  4>  and  0  may  be 
considered  identical,  as  shown  in  the  following  sketch: 


By  reference  to  the  sketch  it  is  seen  that 

d A  2  jt  R*  sin  5 
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Integrating  over  the  front  surface  of  the  sphere,  the  total  number 
of  molecules  directed  upstream  within  do  is 


V-  (Scos0)[l  +  erf  (S 


i  B) 


J  ►  sin  ( 


i  0AB 


«aiR!H| 


—  yjn  S3 
3 


yfv  S5 


erf  S 


(46) 


The  number  of  collisions  per  unit  volume  per  unit  time  at  a  distance 
f  from  the  sphere  is  the  number  of  collisions  per  unit  time  within  do 
from  S.  to  S.  +  d  Si  divided  by  the  volume  of  the  space  within  do  from  i  to 
&  +  d  i  .  This  volume  is  ( i  +  R)Jdcudi  . 


The  expression  for  the  number  of  collisions  between  free-stream  and 
reemitted  molecules  per  unit  volume  per  unit  time  may  be  written  as 


n 


H2Pi 

V  tr  { t  +  r  )a 


—  \f  n  S 


\]  n  S 


erf 


S 


e  — i/A.w 
Aw 
(47) 


Since  the  number  of  free-stream  molecules  and  the  number  of  wall  mole¬ 
cules  colliding  are  equal,  the  number  per  unit  time  of  free-stream  mole¬ 
cules  prevented  from  impact  upon  the  sphere  can  be  determined  by 
integrating  the  above  expression  from  ■  to  the  surface,  4  «  0,  i.  e.  , 


U,*,  -  ni)  rf  d  4 


(48) 


Integrating  Eq.  (48),  using  Eq.  (47)  for  *?,  gives 


( n„ 


erf  S 

3  3 


(49) 


R  /  R  \  R/Xw  f  „./  R  \~l]  f‘2  kT|  \  m 

\k~  vrj; e  r  ,jJu  aJjji  *-)  w 


where  -  Ri  (-R/Aw)  is  the  exponential  integral  (Ref.  36). 

Since  it  is  assumed  that  the  probability  of  a  free-stream  molecule 
being  in  a  collision  is  independent  of  the  thermal  random  velocity  distri¬ 
bution,  then  Ti  is  equal  to  ,  Using  this  equality  and  the  expression  for 
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S,  Eq.  (49)  becomes 

/i  .-s’  VZ  .  V”  erf  s) 


-  ^  .  /  j _ ;  s  t  «  s  .  ±a_  .  ii_ 

n  '  \fi  s*  fiS*  3S  3  3 


From  Eqs.  (45)  and  (50)  it  is  found  that 
cD  *i_  i 


where 


'  -  D 


f(S) 


SM  f-B_  -  ( -JlV  Fa  (- 

nr  \k  \AW/  L  l  OJ 


1  e“S  e  S  1  , —  \fir  ,  c 

- -  +  -  +  —  y  n  +  -  erf  b 

6  S*  fi  $*  3  S  3  3 


And,  from  Ref.  35,  CDfra  is  given  by 

s*  +  /4S4  4  4Sa  -  i 

V  2  S4 


2  Sa  +  l  -  S* 

“*  0  f  n  ”  / 


vTsJ 


2  \T*  rnr 


)  erf  S  +  1±L 

t  3S  \  T 


(50) 


(51) 


(52) 


(53) 


Equations  (51),  (52),  and  (53)  can  be  used  to  obtain  a  drag  coefficient 
of  a  sphere  in  the  near-free-molecular  regime  where  d  «  Aw.  If  the  mean 
free  path  of  the  reemitted  molecule  is  on  the  order  of  the  sphere  diameter, 
then  the  assumption  that  the  sphere  may  be  treated  as  a  point  source  is 
not  justified.  In  the  remainder  of  this  section  an  analysis  is  made  wherein 
the  finite  size  of  the  sphere  is  considered.  It  is  assumed  in  the  following 
analysis  that  the  Mach  number  is  large.  Other  assumptions  are  the  same 
as  those  in  the  preceding  analysis. 

The  number  of  collisions  per  unit  volume  per  unit  time  at  a  point  in 
space  involving  molecules  which  emanate  from  a  surface  element  dA  is 

t/7 

a,.  „-tM-  a  A  (54) 

2„,/>  A.  «’ 

where  nw  is  given  by  Eq.  (17).  For  large  values  of  S  (large  Mach  number), 
nw  becomes 

"w  -  ni  (2  \in  S  cos  0)  (55) 
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In  order  to  get  the  total  number  of  collisions  per  unit  volume  per 
unit  time  at  a  point  P  ,  Eq.  (54)  must  be  integrated  over  the  surface 
elements  of  the  body  that  "see”  the  point  in  space. 

It  is  assumed  that  the  sphere  is  in  a  uniform  free-molecule  flow  com¬ 
posed  of  the  incident  molecules.  The  number  density  of  this  incident  flow 
can  be  determined  at  the  stagnation  point  of  the  sphere.  The  following 
sketch  shows  a  point  on  the  centerline  ahead  of  the  stagnation  point  of  a 
sphere  and  related  notation: 


Since  the  flow  is  axially  symmetric,  the  element  of  area  dA  in  Eq.  (54) 
may  be  written  as 


2  n  R2  sin  9  d  6 

The  collision  density  involving  molecules  reemitted  from  this  area  is 

i/a  — 

dr,  «  Qw  /2kT„N  cos  0  R2  sin  6  g  -i/Aw  ^  $  (56) 

'  m  *  Aw  t  2 


This  expression  must  be  integrated  over  values  of  9  from  0  -  0  to  <£  =  */2  . 
The  total  collision  rate  density  at  P  is  therefore 


V 


J*<p  =  7r/i 
0=0 


C.  OS  R*  BID  9  X/Aw 

Aw  i1  6 


(57) 


In  order  to  evaluate  this  integral,  it  is  convenient  to  express  all  variables 
in  terms  of  &  /Aw  .  From  the  law  of  cosines  for  triangles  with  sides  1,  r,  R. 
it  is  seen  that 

i2  =  r2  +  R2  -  2r  R  cos  9  (58) 
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and 

rJ  =  P2  --  Ra  —  2  £  R  f  —  con  <b) 

Differentiating  Eq.  (58)  gives 

siQ  0  AB  •»  (i/r  R )  d-f 

Let 


and 


It  is  observed  that 


]LJ  Aw  =  A  when  <f>  =  ff/2 

4/Aw  -  B  when  0  «  0 


(59) 

(60) 


(61) 


and 

B  =  r/Xw  -  RAw  (62) 

Substituting  Eqs.  (55),  (58),  (59),  (60),  (61),  and  (62)  into  Eq.  (57)  gives 
the  result. 


■/:  {*(fcX=?)7[fci  i(rj]}[(ti'  -fei  -(£)] 


(63) 


[(*)■  *  (£)'  -  (d]  ‘  (t) 

Equation  (63)  may  be  written  as 

•-  fM  -  ■  (fc)'(t)'  •  (d  -  ■  O) 

“(day.  ai 


V 


Evaluating  the  integral  in  Eq.  (64),  one  finds 


Dj  L'w 


(Ba  +  2B  *  2)  -e  ' 


(A1  +  2.4  +  2)  -  2  (-L\  (65) 

■*-(*)' oa  i 


Equation  (65)  is  presented  graphically  in  Fig.  (12)  where  v  Aw/(m  U^)  ig 
plotted  as  a  function  of  r/Aw  -  R/Aw  for  various  values  of  R/Aw  • 
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Using  Eqs.  (45)  and  (48),  it  is  found  that 


i  r  'aD/ 

rDfm  I  *  __L_  J  „  d  i 

The  integral  in  Eq.  (66)  was  evaluated  numerically. 

Application  of  Eqs.  (66),  (25),  and  (53)  to  calculate  the  drag  coeffi¬ 
cients  corresponding  to  conditions  in  the  LDH  Tunnel  shows  that  Cd  is 
underestimated,  though  not  au  badly  as  when  earlier  theories  are  used. 
However,  it  is  important  to  realize  that  the  greater  part  of  the  experi¬ 
mental  data  lie  outside  the  first-collision  regime.  Comparison  of 
Eq.  (66)  with  the  earlier  theories  represented  by  Eqs.  (29),  (30),  (33), 
(34),  and  (37)  reveals  that  all  are  in  close  agreement  concerning  the 
initial  departure  of  Cd  from  the  free-molecular  limit.  When  lower 
Knudsen  numbers  are  considered,  the  theories  diverge  most  rapidly, 
with  Eq.  (66)  showing  closer  agreement  with  the  data  as  well  as  a  form 
which  is  more  similar  to  the  trend  of  the  data  toward  the  limit  corre¬ 
sponding  to  near-inviscid  flow. 

A  degree  of  uncertainty  exists  because  of  unresolved  questions, 
such  as  those  pertaining  to  the  free-stream  intermolecular  force  field, 
degrees  of  accommodation,  and  manner  of  reflection,  plus  the  neglect 
of  scattered  molecules  and  collisions  subsequent  to  the  first.  For 
example,  the  assumption  of  specular  reemission  results  in  the  present 
theory  being  in  closer  agreement  with  experiment,  as  also  does  the 
assumption  of  less  than  complete  accommodation. 

Following  from  the  observation  of  the  satisfactory  form  of  Eq.  (66) 
and  in  view  of  the  uncertainties  mentioned  above,  the  authors  chose  the 
factor  of  proportionality  in  Eq.  (27)  to  achieve  a  matching  of  Eq.  (66) 
and  the  experimental,  hypersonic,  cold- wall  data.  This  factor  was  found 
to  be  2.  446,  i.  e. , 


Knw  (mod.)  «  2.446  j-- (67) 

The  present  and  all  the  other  theories  reviewed  herein  are  compared 
with  the  new  data  in  Fig.  13  where  Knw  (mod.  )  is  used  as  defined  by 
Eq.  (67).  Use  of  the  factor  2.  446  improves  agreement  of  all  theories 
with  the  data,  but  most  of  the  earlier  theories  still  indicate  much  too 
rapid  a  decrease  of  drag  as  Knudsen  number  decreases.  The  theory 
represented  by  Eq.  (37),  it  will  be  recalled,  already  includes  a  free  con¬ 
stant  chosen  to  best  fit  the  present  data,  but  it  underestimates  the  extent 
of  the  transition  between  near-inviscid  and  near-free-molecular  flows. 
Since  the  factor  introduced  in  Eq.  (67)  and  Fig.  13  probably  may  be 
regarded  as  a  constant  for  hypersonic,  cold-wall  conditions,  its  use  is 
thought  to  be  justified. 
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As  an  aid  in  the  further  examination  of  trie  validity  of  Eq.  (66),  it  is 
compared  with  experimental  data  corresponding  to  2  <  <  li  in  Fig.  14. 

There  it  may  be  seen  that  all  data  are  reasonably  well  correlated  by  the 
chosen  parameters  at  high  Knudsen  numbers.  *  Also,  although  points 
corresponding  to  low  Mach  numbers  are  not  plentiful  at  the  higher  Knudsen 
numbers,  it  may  be  seen  that  Eqs.  (66)  and  (67)  with  the  factor  2.  446 
treated  as  a  constant  may  be  a  satisfactory  approximation  at  the  higher 
Knudsen  numbers,  even  for  Mach  number  as  low  as  two. 


4.0  DISCUSSION  AND  CONCLUSIONS 


First  of  all,  it  would  appear  that  a  semiempirical  approximation  to 
the  drag  coefficient  is  capable  of  closely  matching  the  experimentally 
determined  data  to  nearly  the  free- molecular  flow  regime.  The  provi¬ 
sional  values  of  the  coefficients  in  the  approximation, 

Cd  »  CDi  +  Kt/^/  Re,  +  K,/Re, 

appear  to  be  functions  of  Mach  number  and  heat  transfer  to  the  spheres 
in  perfect  gases.  At  Mach  numbers  above  six,  the  heat-transfer  condi¬ 
tion  becomes  more  important  than  Mach  number.  Values  given  for  Kt 
and  K,  in  the  hypersonic,  cold-wall  condition  are  believed  to  be  more 
accurate  than  those  for  low  Mach  numbers,  largely  because  of  the  stronger 
possibility  of  model  support  and  base  pressure  influences  in  the  latter 
case.  It  is  shown  that  K*  is  negative. 

The  new  first-collision  type  of  analysis  presented  herein  is  shown 
to  produce  an  equation  for  sphere  drag  coefficient  which  possesses 
several  attractive  features.  It  is  not  limited  to  very  high  Mach  numbers 
and  cold  walls,  and  it  yields  a  drag  coefficient  at  the  limit,  Kn  -»0,  which 
is  on  the  order  of  the  known,  continuum,  near- inviscid  value.  Use  of  an 
experimentally  determined  factor,  which  probably  may  be  regarded  as  a 
constant  for  high  Mach  numbers  and  cold  walls,  causes  the  present  theory 
to  fit  experimental  data  satisfactorily  at  Reynolds  numbers  far  beyond  the 
first-collision  regime.  Equation  (66)  with  the  modification  of  Eq.  (67)  is 
shown  in  Fig.  15  to  effect  the  correspondence  of  the  continuum  to  the  non¬ 
continuum  analytical  results,  thus  offering  a  means  for  determining  drag 
of  spheres  under  the  conditions  specified  and  throughout  the  entire  range 
of  Reynolds  or  Knudsen  numbers  encountered. 


*The  data  would  approach  different  limits  as  Knudsen  numbers 
approached  zero,  since  Cnj/Cnf[D  *  const,  for  2  <  <  11. 
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TABLE  1 

CONDITIONS  AND  RESULTS  OF  THE  LDH  WIND  TUNNEL  FREE-FALL  TESTS 


h'yy  /  hg 

M. 

Diameter, 

in. 

Re*, 

Re  2 

CDfm 

CD 

CD<'CDfm 

0.  076 

2.00 

10.  6 

0.  0312 

9.  40 

1.08 

2.  21 

1.  94 

0.  877 

2.00 

10.  6 

9.  46 

1. 09 

2.  02 

0.  912 

2.  00 

10.  6 

0.  0392 

11.  95 

1. 38 

1.  99 

0.  899 

2.00 

10.  6 

11.  86 

1.37 

1.  94 

0.  874 

2.  00 

10.  6 

11.  86 

1.  37 

1.88 

0.  849 

1.  98 

10.  5 

0.  0469 

14.  46 

1.  69 

1.  90 

0.  859 

2.  00 

10.  6 

14.  19 

1. 64 

1.  94 

0.  874 

2.  00 

10.  6 

0.  0625 

18.  92 

2.  19 

1.  87 

0.  846 

2.  05 

10.  7 

0.  0781 

23.  07 

2.  64 

1.  76 

0.  794 

2.  01 

10.  6 

23.  49 

2.  71 

1. 78 

0.  805 

0.  086 

2.  29 

10.  7 

0.0312 

10.  74 

1.  18 

2.22 

2.  03 

0.  914 

2.  25 

10.  6 

11.  01 

1. 23 

1.  95 

0.  878 

2.  26 

10.6 

10.  95 

1.  22 

1.  94 

.  0.874 

2.  28 

10.  7 

0. 0392 

13.54 

1.  49 

1.  92 

0.  862 

2.  26 

10.  6 

13.  74 

1.52 

1.  98 

0.  891 

2.  28 

10.  7 

0. 0469 

16.  20 

1.  78 

1.93 

0.  868 

2.  24 

10.6 

16.  60 

1.  86 

1.  94 

0.  870 

2.  26 

10.  7 

0. 0625 

21.  81 

2.  41 

1.  85 

0.  834 

2.  28 

10.  7 

21. 59 

2.  37 

1.  09 

0.  896 

2.  24 

10.  6 

22.  13 

2.  40 

1.  88 

0.  847 

2.  26 

10.  7 

21.  81 

2.  38 

1.  86 

0.  838 

2.  24 

10.  6 

22.  13 

2,48 

1.  82 

0.  820 

2.  29 

10.  7 

0.  0781 

26.  85 

2.94 

1.  68 

0.  756 

2.  24 

10.  6 

27.  95 

3.  14 

1.  69 

0.  761 

2.  28 

10.  8 

0. 0938 

31.  89 

3,47 

1.  65 

0.  741 

2.  30 

10.  7 

32.  06 

3.50 

1.  68 

0.  753 

2.  30 

10.  7 

32.  06 

3,50 

1.  69 

0.  762 

2.  26 

10.  7 

32.  72 

3.62 

1.69 

0.  758 

2.  24 

10.  6 

33.  37 

3.  74 

1. 62 

0.  728 

2.  24 

10.6 

33.  55 

3.  77 

1.  63 

0.  732 

2.  25 

10.6 

33.  03 

3.68 

1.  72 

0.  774 

2.  24 

10.  6 

33.  55 

3.  77 

1.  61 

0.  723 

31 


AEDC-TDR-62-205 


TABLE 

1  (Concluded) 

h-yv  /  ho 

Tw/T. 

M. 

Diameter, 

in. 

Re. 

Fte2 

CD 

Cd/C  Df, 

0.  086 

2.30 

10.  7 

0.  1250 

42.  75 

4.  66 

2.  22 

1.  60 

0.  717 

2.  30 

10.  7 

42.  75 

4.66 

1.  68 

0.  754 

2.25 

10.  6 

44.  04 

4.  81 

1.  56 

0.  703 

2.  29 

10.  7 

0.  1562 

53.  72 

5.88 

1.55 

0.  696 

2.  25 

10.  6 

55.05 

6.  14 

1.56 

0.  700 

2.  24 

10.  6 

55.61 

6.23 

1.  61 

0.  724 

2.  29 

10.  7 

0.  1875 

64.45 

7.05 

;.  55 

0.  606 

2.  31 

10.  8 

0.  2188 

74.07 

8.03 

1.  50 

0.674 

2.  29 

10.  7 

0.2500 

53.70 

9.  40 

1.48 

0.  667 

2.  26 

10.  7 

0.2812 

98.  16 

10.86 

1.  44 

0.  646 

2.  30 

10.  7 

0. 3125 

106.87 

11.66 

1.  42 

0.  638 

2.  30 

10.  7 

0.3438 

117.56 

12.  82 

1.  42 

0.  638 

2.  28 

10.  7 

0.3750 

129.57 

14.  21 

1.  36 

0.  609 

0.  121 

3.03 

10.  5 

0.0312 

13.65 

1.  38 

2.  26 

2.  07 

0.  918 

3.  03 

10.  5 

13.65 

1.  38 

2.  00 

0.  886 

3.  06 

10.  6 

13.  42 

1.  35 

1.  89 

0.  839 

3.  03 
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Fig.  1  The  Low-Density,  Hypervelocity  Wind  Tunnel 
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Fig-  2  Typical  Trajectory  of  Sphere  during  Free  Flight  through  Test  Section 
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Free-molecular  limit  for  diffuse,  fully  accommodated  re-emission: 
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Fig.  5  Drag  of  Spheres  in  Hypersonic  Flow  with  <  h0:  New  Dotni 
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Fig.  8  Provisional  Values  of  K] 
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Fig.  9  Provisional  Values  of  K2 
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Fig.  10  Nomenc lature  Used  lor  Incident  and  Reemitted  Molecules 
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Fig.  13  Comparison  of  Theories  and  LDH  Wind  Tunnel  Data  for  the  Hypersonic,  Cold-Wall  Condition 


